A Green's function approach to transmission of massless Dirac 
fermions in graphene through an array of random scatterers 

Neetu Agrawal (Garg)-*^, Sankalpa Ghosh^ and Manish Sharma-^ 
^Centre for Applied Research in Electronics^ 
Indian Institute of Technology Delhi, New Delhi-110016, India and 
'^Department of Physics, Indian Institute of Technology Delhi, New Delhi-110016, India 

(Dated: March 6, 2013) 

Abstract 

We consider the transmission of massless Dirac fermions through an array of short range scat- 
terers which are modehed as randomly positioned 5- function like potentials along the x-axis. We 
particularly discuss the interplay between disorder-induced localization that is the hallmark of a 
non-relativistic system and two important properties of such massless Dirac fermions, namely, com- 
plete transmission at normal incidence and periodic dependence of transmission coefficient on the 
strength of the barrier that leads to a periodic resonant transmission. This leads to two different 
types of conductance behavior as a function of the system size at the resonant and the off-resonance 
strengths of the delta function potential. We explain this behavior of the conductance in terms 
of the transmission through a pair of such barriers using a Green's function based approach. The 
method helps to understand such disordered transport in terms of well known optical phenomena 
such as Fabry Perot resonances. 

PACS numbers: 72.80.Vp,71.23.An,73.23.Ad, 72.10-d 
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I. INTRODUCTION 

The remarkable transport properties of graphene are primarily due to the fact 

that under ambient conditions the charge carriers are massless Dirac fermions with definite 
chirahty. Such electrons get differently scattered by a potential barrier [3] as compared 
to non-relativistic electrons in other conventional semiconductor or metal . This led to a 
number of transport anomaly in the graphene. The expression of conductivity of ballistic 
graphene is remarkably diflFerent from that of a non-relativistic electron and the mini- 

mum conductivity is given ^ for undoped graphene. The pref actor 4 is due to contribution 
from two sublattice degrees of freedom and two valley degrees of freedom. 

This prediction was verified experimentally [6|,17l]. More remarkably, addition of impurities 
to pristine monolayer graphene leads to the conductivity enhancement at the Dirac point, 
where as addition of such impurities away from the Dirac point leads to a supression of 



conductance 



10l |. A general theory to understand the transport of Dirac fermions in 



presence of such isolated impurities was also developed [11] A more complete theory that 
takes into account the effect of disorder as well as interaction effect was also constructed. 
I12I that gives better agreement with the transport measurement. Such impurity scatterers 
could be vaccancies, adsorbed atoms, molecules, or impurity clusters, {xstll^, or hydrogen 
atoms controUably added to the surface [15] or metallic islands deposited on graphene surface 
{xot-The difference between the infiuence of short range scatterers and long range Coulomb 
impurities on the transport of massless Dirac fermions was also studied |l7|. Theoretical 
progress was also made to understand the nature of transport in presence of correlated 
disorder [isl . 

For a better theoretical understanding of these transport properties of massless Dirac 
fermions in graphene, it must be compared with the conventional charge transport in dis- 
ordered condensed matter system which is principally understood in the frame work of the 
pioneering work of Anderson[19] and subsequent developements [20] in this direction. Par- 
ticularly important in the context of graphene which is a two-dimensional atomic crystal 
the prediction of scaling theory of localization [21]. The scaling theory predicts that below 
two dimension any amount of disorder can localize all the states. One of the implications 
of this scaling theory is that conductace G approaches zero as the sample size L goes to 



infnity for a disordered one dimensional system and such decay is exponential in nature 



22|. 
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These predictions are based on the fact that charge carriers are non-relativistic in nature 
and obeys Schrodinger equation. Thus it is extremely important to study the revision of 
above weU estabhshed properties of non relativistic electrons for the case of massless chiral 
dirac fermions that dominates the transport properties of graphene. 

The transmission of such massless dirac fermion through one dimensional potential barri- 
ers demonstrates two fundamentally different behavior as compared to the similar transmis- 
sion of non relativistic electrons First of all, they Klein tunnel through such barrier which 
implies full transmission at normal incidence. Secondly the transmission prbability periodi- 
cally oscillates with the varying strength of the barrier. This particular property which was 
already implicit in the transmission expression given in jsl was more clearly demonstrated 
in the ref. [8] by considering transmission through short range scatterers approximated as 
delta function potential. 

The above mentioned exotic properties of transport electrons led to a large body of work 
devoted to the study of ballistic progation of two dimensional massless Dirac fermions in 



graphene t 



tentials 



irough various one dimensional superlattices made out of scalar and vector po- 
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28| . Experimentally also certain superlattice structures were imposed on mono- 
layer graphene j^]. Relatively much less attention has been paid to similar studies 
through one dimensional array of disordered or impurity potential. In this paper we study 
the transmission of such massless dirac fermions through a one dimensional arrangements 
of short range impurities modelled as delta function potetials on random locations using 
a Green's function based technique. We obtain an analytical expression for the transmis- 
sion and zero temperature conductance through such one dimensional arrangement of delta 
function scatterers. 

Using this theoretical framework we show that the transmission and conductance in the 
presence of array of disorder potentials can be understood in terms of the transmission 
through a double barrier structure consisting of such delta function like potentials. One of 
our main results that the conductance properties through such barriers can have two type 
of behaviors as a function of the system size. If the strength of the each delta function 
barrier is close to its resonant values js^ then there is perfect transmission through each 
barrier and the conductance does not change significantly as the system size increases. On 
the otherhand if the strength of such delta function scatterers is well oS from the resonant 
value, the conductance shows a algebraic decay as a function of the size of the system with 
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an exponent which is evaluated to be close to 0.5. As our results shows this happens over 
a wide range of the strength as well as the mean separation between two such successive 
scatterers. The Green's function method that we use is quite general and can also be used to 
reporduce well known results (25 . 26| for ordered superlattice such as Kronig Penny systems 
and their different variants [32l|. 

The paper has been organized in the following manner. In various subsections of section 
im we develope the theoretical framework of the Green's function method and derived the 
expression of transmittance and conductance. In Section [TTIl we first explain the single and 
two barrier transmission and explain how the barrier transmission can be understood in 
terms of successive two barrier transmission. We then present the results for transmission 
through N such barriers, first when the positions of such barrier is randomly located, but 
all having equal strength, and then with both position as well as strength randomized. We 
conclude by pointing out a possible generalization of our technique and implications of our 
result in Section I 



II. GREEN'S FUNCTION AND THE TRANSMISSION 

In this section we first obtain the free particle Green's function for massles Dirac fermions 
and determine the solution in presence of short range scatterers modelled as an array of the 
delta function potentials. We then analyse how these solutions can be understood in terms of 
multiple scattering processes takes place between two such barriers. We shall finally obtain 
the expression for the transmittance through A-such barriers with arbitrary position and 
strength. 



A. Wave functions of massless dirac fermions in terms of free particle Green's 
function 

The charge carrier in Graphene under ambient condition behave like two dimensional 
masssless Dirac fermions [3]. For such charge carriers with energy the stationary solutions 
are obtained from the following Dirac- Weyl equations 

— ihvF^ ' V'0(x, y) + Vlip{x^ y) = Eil){x^ y) (1) 
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Here ^ 10^ms~^ is the Fermi velocity and a = (JxX+cjyy. Here y) is a two component 
pseudo-spinor where the pseudospin refers to the sublattice degrees of freedom. We take 
V{x^y) = V{x) such that translational invariance along y-directions leads to ip{x^y) = 



^ikyv ^ Rescaling lengths and energy like quantities by Fermi wavelength Ai? = 



27r 



and hvp/ls) we introduce dimensionless variables x = x/\f-> = kx^yXp^ e = EXf/Hvf^ 
and V = VXp/hvp. Substituting the ^-translationally invariant solutions in Eq. ([T]) and 
multiplying by a^;, in terms of the dimensionless variables mentioned, we get the effective 
one dimensional equation 

d _ I —iky e \ _ _ 
^^^(^) + _ _ ^(^) = v{x)axilj{x) (2) 

We are interested in seeking the nature of transmission of such massless Dirac fermions 
through a potential of the form 

N 

^(^) = - xi) (3) 

1=1 

The above potential is a series of delta functions that are randomly positioned within a length 
(say) L, such that the end points are held fixed at Xi = and Xn = L. A delta function is 
placed at each of the edges Xi and Xn while the number in between the edges can be varied. 
Additionally, the strength of each of the delta function can also be varied. Thus we are 
studying the transmission of two dimensional Dirac fermions through an one dimensional 
potential that mimics the eflFect of a set of random short range scatterers on a line. If the 
range of the potential of isolated impurities is much smaller than the fermi wavelength of 
electrons ( which is theoretically infinity at the Dirac point), but larger than the carbon- 
carbon bond length in monolayer Graphene, the scattering potential of an isolated impurity 
can be modelled as a delta function. 

Here we describe the main steps of the Green's function based method that we adopted 
to calculate such transmission. The method particularly takes care of the massless ultra 
relativist ic nature of such fermions. 

The eq.([2]) is an inhomogenous differential equation whose full solutions can be written 
as a sum of the complementary fnction and particular integral of which the later can be 
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expressed in terms of the Green's function. The resulting solution can be written as 

/oo 
dx'G{x,x'Ya^v{x')ij{x'), (4) 
-OO 

where the Green's function xY is a 2 x 2 matrix given as 

G{x,x%^-, = e'~'^^^-^'^G'' (5) 
The pseudopsin dependent part of the Green's function is given by 




G'' = 7r7zn\ ,\ (6) 

Similarly 



1 



G{x,x%^,, = I ' " .. I (7) 




In the above solutions 74 = ^41 | , 5 = Si | I , and Bi are constants. 

(j) = tan~^ky/kx with kjp = kl + ky 
Substituting the form of potential given in (|3j) in the solution (|4]) we get the stationary 
solution at any spatial point x as 

^[x) = A/^^^ + Be-'~^^^ + XiG{x, xiya^iPixi) (8) 

1=1 

B. Multiple scattering processing and the solution at N-th barrier 

The above solution ([8j) implies that we need to determine ip{xi) at the locations for each 
delta function, namely at / = 1 to N to get the solutions at any spatial point x. This is due 
to the fact that they are the only scattering centers that can change the incident amplitude. 
It may be pointed out that above method of calculating the amplitude through the solution 
dlj) can be implemented for any other form of potential in one and higher dimension as well. 
However in the presence of delta function like scattering centers the integral collapses into 
analytically tractable series of the series given in (|8j). Finally after having the solution at each 
point of space we can calculate the ratio of transmitted current density and incident current 
density to yield the transmission coefficient through such series of scatterers. However before 
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FIG. 1: Schematic display of the matching condition in presence of delta function barrier. 

proceeding to determine the expression for the transmission coefficient using this method 
we shall discuss briefly some issues that are relevant to the determination of the value of 
i/j{xi) at the location of the dirac delta function. 

For a non-relativistic spinless particle obeying one dimensional Schrodinger like equation, 
the solution at the location at the delta function is obtained by taking the average of the 
left and right hand limit, namely 

^(0) = l[^(0+) + ^(0-)] 



This situation is schematically pointed out in Fig. [H However as pointed out in a number of 



works [25|, 



- 35| this procedure cannot be used in determining the solution at the location 
of the delta function for the corresponding relativistic problem obeying the one dimensional 
dirac equation which is a flrst order diflFerential equation. The issue is resolved by matching 
the two component spinorial wavefunction on the right and left side of a flnite width barrier 
and then flnding out of the form of the transfer matrix by allowing such fnite width barrier 
to approach delta function limit. 
Such matching condition gives is 



>2(0+) 



>2(o-) 



(9) 





( COS A 


—i sin A \ 










I — isin A 


cos A / 
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We note that S is a periodic function of A and that S = ±7 for A = nvr. This the special 
situation for which the barriers behave as if they are transparent to the incident charge 
carriers. 

Using the above form of the matching matrix, the wavefunction at the position of the 
delta function is expressed as the following linear combination: 

ilj{xi+) + ^p{xi-) = A cot ^^p{xi) (10) 

A 7^ 2n7r, n = 0,1,2,3--- js^. The corresponding method the case of A = nvr where the 
matching matrix is a unity matrix and leads to full transmission through the delta function 
scatterers via the matching condition (|9j) and will be described later. The above method will 
now be implemented first to determine to evaluate for all 1 values, namely for / = 1 

toA^. 

To calculate the transmission through such delta function barriers we need to evaluate 
the wavefunction at the last such barrier which is the A^-th one, namely iP{xn) • According 
to eq. [8] this can be written in terms of i/j{xn±). Now, 

N 

^P{xn+) = Ae'~^-^^ + Be-'~^^^^ + A^G(x^+, xiYa^^ixi) 

1=1 

= {A + S_{N))e'~^-^^ + 5e-^^-^^ (11) 
Here we used the scattering matrix S-{N) defined as 

S-{N) = exp{-ik:,XN)Y^XiG{xN+,xiya:,ilj{xi) (12) 

1=1 

which gives us the modification of the amplitude in the forward propagating wave A due to 
scattering by all the A^ barriers. 

Similarly the amplitude at the immediate left of the A^-th barrier 

N 

^{xn-) = ^e^^-^^ + Se-^^-^^ + \iG{xn-, xiYa^^ixi) (13) 

1=1 

Using the above expressions and the relation ( |TOl) a straightforward calculation can de- 
termine 



^{xn) = exp{ia • x^) {A + S^N)) e^^-^^ + Se"^^-^^ (14) 
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The above expression can be understood as follows. The first term ^j^^ which provides 
the envelope of the propagating wave through the delta function barrier is due to the Fraun- 
hofer diffraction of the wave function by the delta function barrier. The second term is due 
to that because of the linear dispersion of the massless Dirac fermions, the delta function 
potential barrier directly changes the phase of the pseudospinor by changing the x compo- 
nent of the wave vector. This fact is expressed by the explicit presence of the pseudopsin 
rotation operator, where the angle through which the pseudospin rotation takes place is 
directly proportional to the strength of the delta function potential. The third term is the 
modification of the free particle wave function due to the multiple scattering from the 
delta function barriers. It may be noted because of the presence of S-{N)^ the scattering 
matrix the determination of iP{xn) requires the determination of for / = 1 to A^. The 

expression for ip{xi) can again be obtained in the same way as in the case oi iP{xn)^ yielding 

ij{xi) = ^ exp(^<T ■ x^) {pi[A - S.iN)y~'^^^ + piBe-'^^'-^] (15) 

2 

Since the right hand side of the above equation also contains to determine the the 

scattering matrix S-{N)^ it has to be solved self consistently to get the solution ^l){xi) 
The term p which is given as 

J2-1 J3-1 J2-1 

Pi = I+ Vjii + J2 Yl ^ni^nn + Y Y ^Jii^J3j2Vj2ji + ' ' ' (16) 

j=l+l 32=1+2 ji=l+l j3=M 32=1+2 h=l+l 

Where 

^ji = A^-^-^exp(zfc^(% - xi))[G{xj,xi)l^^^^ - G{xi,Xj)l^^^^]a^exip{ia • x^) (17) 

2 

and fjji is just the complex conjugate of the r]ji. 

Thus rjji transfer the amplitude from j the delta function to the / delta function due 
to the scattering at the j-th delta function. The series pi indicates how the amplitude is 
transferred to the l-th barrier from all the barriers that succeeds it as one goes from left to 
write through the multiple scattering processes of various order. 

The first term is an identity matrix and indicate the zeroeth order process, namely when 
I = the final barrier as was the case for the expression given in f fMl ). The second term 
refers through a single scattering process by which the amplitude is transferred from all 
the barriers right to the l-th barrier to the l-th barrier through a single scattering. The 
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third term refers to the processes where such amphtude transfer takes place through two 
successive scatterings and hence a second order process. 

The scattering matrix S-{N) can be now be self consistently determined by substituting 
eq. [15] into eq. [12] which can be written in terms of pre determined quantities such that 

-1 



S-iN) = 



N 



1=1 



N 

1=1 



(18) 



where 



Xi 



V- 



G ax exp[ia • x — 



(19) 



where is the spin compoent for the Green's function . 



C. Calculation of Transmission coefficient 



In our present theoretical framework we are assuming the massless dirac fermion to 
incident from the left side of the first barrier and after A^-the barrier it will again freely 
propagate towards the right. The structure of S-{N) and pi implies this assumption. To 
calculate the transmission coefficient in the presence of (say) N number of barriers we need 
to evaluate the wavefunction for x < Xi and for x > Xa^. For x > Xa^, all the scattering 
centres which lie behind x^ contribute. Also since there is no wave coming from the right 
side we set B = 0. Hence 



ij{x) = [A + S-N]e'~^^^ x> xn 

The transmittance T can be obtained as the ratio of the transmitted probability density 
and the incident probability density is obtained as: 

T = \tf = [A + S_N]-^[A + S_N]/A^A 



Setting the condition B = in the expression in the expression f [T8l ) of the scattering 
matrix S-{N) 



A + S^N) = 



1 



N 



det 



I + 22i=i XiPi 



-cof 



A 
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Using the relation [36] 



N 



CO 



^factor I + ^XiPi A = A 



1=1 



we obtain 



2 



1 



(20) 



T 



det I + Y.i=iXiPi 



It may be mentioned for A = nvr, T = 1 since the matching matrix is an identity matrix. 
The above expression provides us for a given angle of incidence 0, the transmittance through 
a randomly positioned delta function like barriers whose strength can also vary from 
point to point. The expression for the two terminal conductance can now be obtained 
by suitably integrating the above expression for all possible angle of incidenc using the 
expression G = Gq JJ^^^T(A, 0) co^(j)d(j). Here, Go = ^EpLyC^ / {vph^) and Ly is the width of 
the system. 

For the current problem the expression for two terminal conductance in the presence of 
random delta function like impurities on a line will be 



It maye be mentioned that the same problem can also be solved with the transfer matrix 
method. But the interpretation of the finally transferred amplitude in terms of a multiple 
scattering process at various orders comes out explitly in the Green's function framework. 
This provides a way to understand scattering through random barriers with the help of 
scattering through two barriers. A connection between the Green's function method and the 
transfer matrix method for non-relativistic electrons is provided in ref . [36] . The method is 
also similar in spirit to the invariant embedding approach used in radiative transfer processes 



In the subsequent section we use this expression to evaluate the transmittance and con- 
ductance under various conditions and analyze these results. 
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(21) 



37]. 



III. 



RESULTS AND DISCUSSIONS 



In this section we shall shall use the expression (1^ and to calculate the transmis- 
sion and resulting conductances through such various combinations of the delta function like 
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barriers. We shall first analyze the case of the transmission through a single delta function 
like barrier and two delta function like barrier to understand some peculiarity of one dimen- 
sional transmission of two-dimensional massless Dirac fermions in Graphene through such 
single and double barrier structure as compared to a corresponding non relativistic prob- 
lem. We explain how these peculiarities are expected to change the localization properties 
of such two dimesional massless Dirac fermions. Subsequently we shall present our result 
for transmittance and conductance through a random array of such delta function barriers 
over a wide range of strength as well as delsity of such barriers and discuss the nature of 
localization of such massless Dirac fermions in presence of short range scatterers. 



A. Single delta function barrier 

The transmittance of two dimensional massless Dirac fermion of energy E [s^ through a 
one dimensional potential barrier of height Vq and width D such that |Vo| >> is given by 
the well known expression 



cos^ 



T - (22) 

1 — cos^{qxD) sin^ 



where Qx = y^^r — ^ and ky = The above expression has two important diflFerence 

with the corresponding non-relativistic problem. At normal incidence for = 0, T = 1 due 
to Klein tunneling. Also everytime the condition QxD = nn is satisfied T = 1, showing 
resonant transmission with tt periodicity. 

The delta function limit of such square barrier is taken by allowing Vq ^ oc and D ^ 
, such that A and the resulting transmission can be obatined as 

T 



1 + sin^ A tan^ (j) 

The above expression can be directly obtained from the expression (l20l) after setting in the 
equation [3] 

v{x) = XS{x) (23) 

A comparison with the corresponding results for a non-relativistic problem shows that 
in the case of the later the transmittance decays with the increasing strength of the delta 
function barrier. 
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B. Two delta function barriers 



Next we study the case of transmission of such massless dirac fermions through two such 
delta function barriers which are separated by a distance d ( in dimensionless form). The 
scattering potential takes the form 

v{x) = XiS{x + ^) + ^2S{x - ^) 

Such a potential structure forms a Fabry perrot cavity for the massless Dirac fermions and 
the transmission through such cavity shows Fabry Perrot resonances [38] . According to ( !20l) , 
for Ai = A2 = A, the transmittance through such double barrier structure is given by 

T = [1 + tan^(l){co8 kxd sin 2X — 2 sin fc^^dsin^ A/ cos(/))^]~^ 

The above expression clearly demonstrates that perfect transmission takes place respectively 
when = 0, A = nvr and cos cj) = tanfca^JtanA. The first two conditions respectively corre- 
spond to the Klein tunneling through such barrier and resonant transmission as the strength 
of barrier is varied. Both features are single barrier transmission feature and reoccurs when 
two such barriers are placed side by side. The third condition refers to the Fabry Perrot 
resonance condition is a purely double barrier feature and occurs due to multiple refiections 
from the two-barrier structure. 

In section ini we developed the expression for transmission through such delta function 
barriers, by writing the scattering matrix S-{N). This in turn is determined by the series 
expression of the pi matrix given in f flGl) whose various terms depicting the different order 
multiple scattering processes. Each such term can be written as a product of ij matrices 
given by the expression f fTTI ) which depicts the transfer of amplitude from one barrier to 
the other. Thus the features through N delta functions barriers, that have the strength, 
but located randomly on a line can be understood through in terms of multiple two barrier 
transmissions. This is one of the main results of the current work. 

The above statement can be well explained using theory of Fabry-Perrot resonances of 
massless Dirac fermions [38] through such double barrier structure. If the strength of the 
delta function barriers are equal and their respective transmission and refiection amplitude 
is given by t and r the multiple scattering processes that involves any two such barriers is a 
linear superposition of terms like tt, ttr^ exp(iC), • • • , exp(i2(n — In all these terms t is 



13 



independent of the separation between two such barriers where as the phase change ( that is 
acquired at each scattering between two barriers is dependent on the separation between such 
barriers. Since for massless Dirac fermions t^(A) = t^(A + nvr), such a hnear superposition 
is also going to be 7r-periodic in the strength of the barrier A. The same argument can be 
extended to the multiple scattering process that involves transfer of amplitude from j-th 
barrier to the i-th barrier through any number of intermediate barriers. Such a transmission 
can be wriiten in terms of the product of the rj matrices defined in f lTTI) each of which is 
going to be periodic in terms of the strength of the barrier. However now these barriers 
being positioned randomly on a line, the phase change induced by each of these pair of 6 
function barriers is random. Thus the total transmission through such barriers will be given 
by a sum of terms whose amplitude is periodic in the strength of the potential A, but phase 
is random. The resulting transmission will therefore show periodicity as a function of the 
barrier strength A, unless the random phase factors will add-up to zero to give a complete 
destructive interference and leads to Anderson localization. The situation is similar to the 
.^ahiy Perrot resonances in disordered one dimensional array of alternating dielectric bi layers 
39I . Such a destructive interference related localization takes place in the related problems 
of non-relativistic electrons governed by Schrodinger equation. However for massless two 
dimensional Dirac fermions such a situation is averted because of Klein tunneling, since 
this ensures that there will be always full transmission at normal incidence. The sharp 
contrast between the transmittance through such one dimensional delta function barriers 
for two dimensional massless dirac fermions and two dimensional non-relativistic electrons 
obeying Schrodinger equation is depicted in Fig. [2l The almost full transmission around 
normal incidence for massless dirac fermions can be contrasted again the exponential decay 
of transmission at any angle for non-relativistic electrons. 

The conductance in presence of such double barrier structures is plotted for such massless 
Dirac fermions in Fig. [3] (b) and compared against the similar double barrier conductance 
for non relativistic particles obeying Schrodinger equation Fig. [3](a). For the massless 
dirac fermions, the transmission shows periodicity as function of the separation between the 
barriers ( Fabry Perrot transmission) as well as the strength of the barrier. The double 
barrier transmission for the non-relativistic electrons that is depicted in the lower figure 
shows periodicity as a function of the separation between the barriers only. There is no 
periodicity as a function of the strength of the barrier for such non relativistic fermions. 
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FIG. 2: (color online) Transmission through disorder for (a) Carriers obeying Dirac hke equation: 
The effect of Klein tunnehng. (b) Carriers obeying Schrodinger equation. The y axis plotted the 
angle of incidence. The x-axis plotted the length of the sample L in the unit of Fermi wavelength. 

since this is attributed to the ultra relativist ic nature of the charge carriers in graphene. 

Next we shall consider the transmission of such massless dirac fermions through a double 
barrier structure when the strength of the barriers are unequal, i. e. Ai 7^ A2. With the 
help of the expresssion f l20l ) we can evaluate the resulting conductance which is plotted in 
Fig. [3] (c) as function of the mean strength of the two barriers along the x~axis as well 
as the difference of strength of the two barriers along the y-axis. Here the amount of 
pseudopsin rotation imparted by the two barriers are of different. Thus the transmission ( 
conductance) resonance occurs when mean Ai + A2 = 2n tt and Ai — A2 = 2n7r, both the 
conditions are satisfied. Also because of the mismatch of the barrier height, the resonace 
peak has a double hump structure which characterizes the resonace due to a double barrier 
structure with a finite difference between the height of these two barriers. It can be now be 
seen since the transmission through many such barriers with random position and strength 
can be thought as a product of transmission through such double barrier structure with 
unequal strength the resonance condition for each such pair will be diffferent from another 
in general and consequently the height of the conductance resonance peak will get reduced 
with increasing fluctuation around the mean strength of such delta function barriers. This 
fact is demonstrated in Fig. [H In Fig. [41(a) we have plotted the dimensionless conductance 

15 




FIG. 3: (color onlinejTwo barrier conductance for non relativistic electrons when the barrier 
strength are equal, (b) Two barrier conductance for graphene charge carriers (b) when the barrier 
strengths are equal (c) Two barrier conductance when the barrier strengths are unequal. 

in the presence of randomly positioned barrier all having the same strength. As one can 
see the periodic occurrence of conductance resonance as function of the strength of the delta 
function barrier. The lower plot of the same figure shows how differential conductance varies 
as a function the strength of the delta function. In Fig. [4](b) we plot how this resonance 
peak changes when apart from the randomness in position we also introduce randomness in 
strength. In the second case we plot the conductance as a function of the mean strength and 
the fiuctuations around this mean. One can see a conductance peak is still observed now 
at the same mean value of A , but with increasing fiuctuations around the mean the height 
of the conductance peak gets reduced. Current experiments can measure the conductance 



and difiFerential conductance in ballistic regimes for graphene based microstructures 



4], 
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FIG. 4: (color online) (a) Conductance as a function of strength of potential in the presence of 
many barriers with position disorder but same potential strength, A = 2. Average separation is 
such that kpL = 5. (b) The change in the shape of conductance resonance as the fluctuation 
around the mean strength of the barrier is allowed to grow. 

Thus our predictions can be directly verified. 

C. Transmission and conductance through N barriers 

1. Barriers with equal strength 

With the above analysis of transmission through a double delta function like structure, we 
shall now analyze the transmission and the resulting conductance through such randomly 
positioned barrier having equal strength A. For a given such strength, we have varried the 
length of the sample L, mentioned in the unit of Fermi wave vector, , keeping the mean 
separation between the disorder /, or the dimensionless quantity kpl constant. Since in all 
our calculations are done for a particular enegy that we call Fermi energy which is in the 
dimensionless unit fci?, different kpl correspond to different mean separation between the 
disorder. The higher value oi kpl thus implies a weakly disordered system where as lower 
value of kpl implies a relatively strongly disordered system {42]. 

With the help of expression f l20l) and f l2T]) we evaluate the dimensionless conductance as 
a function of the sample size, for different strength of the randomly positioned delta 
function potentials for a given kpl and a given incidence energy e. Such results are presented 
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FIG. 5: (color online) Decay of Conductance with increasing sample size, on linear, log-log and 
semilog scale. Average separation between the barriers is such that kpL = 5, e = 1 

in Fig. (Ej) and Fig. (|6j) for two different values of kj^l for a set of disorder or delta function 
barrier strength A. For a given Uength L , G{L) is evaluated after doing ensemble avearge 
of randomly positioned delta function like barriers in that length. 

The plot for dimensionless conductance as a function of L shows that for smaller 
sample size conductance shows fluctuations. This fluctuations are associated with the con- 
ductance fluctuations in mesoscopic samples occurs due to the inhomogeneity in the position 
of the scatterers in such sample and well studied in the literature j43| . We shall not discuss 
this issue further and focus on the behavior behavior of G{L) for larger sample size when 
such fluctuations die down. We found that the L dependence of G{L) can be broadly divided 
in two parts. For the strength of the deta function barriers satisfying resonant transmis- 
sion, namely A = nvr, T = 1 and the conductance remains constant as a function of length. 
Close to this resonant strength, conductance thus show a very slow decay and continuous 
achieve this resonant behavior. This can be straightforwardly understood with the preceed- 
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FIG. 6: (color online) Decay of Conductance with increasing sample size, on linear, log-log and 
semilog scale. Average separation between the barriers is such that kpL = 1, e = 1 

ing discussion of analyzing transmission random delta function like barrier in terms of 
two barrier transmission. Away from this resonance strength, the conductance shows an 
algebraic decay as a function of the sample length L. To extract this algebraic decay for 
different kpl value, we have plotted the log as a function of L as well as logL. Fitting 
these plots we find that for non-resonant strength of the delta function barrier G{L) can be 
well approximated by the expression For the non resonant values of strength of potential we 
obtain: 

= a = 0.46 - 0.56 (24) 
Go L« 

with c being a sample dependent constant. Our results agrees well with the observation in 
ref. flO] where a random-matrix theory based argument also predicts an algebraic decay of 
the conductance with exponent 0.5. 

Finally we plot the conductance on a log-log scale for several values of kpl starting from 
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FIG. 7: (color online) Position as well as potential strength disorder: Linear fit on a log-log 
scale-showing algebraic decay 

the 5 ( weak disorder) to kpl = 0.2 in Fig. [71 The smaller values of kpl correspond to fairly 
strong disorder and one needs to calculate the ensemble averaged tranmsission through a 
very large number of delta function barriers. The fitting of of the plots suggest that the 
algebraic decay again depicts the correct dependence of the conductance on the system size, 
when the later is large. However the conductance fluctuationation persists over a larger 
length scale as the value of kpl is lowered. The much slower decay of the conductance 
near the resonant value the mean impurity strength also get more and more supressed as 
suggested by Fig. [H (b) with increasing fluctuations around the mean strength. 

IV. CONCLUSION 

We conclude by summarizing our main findings. First the Green's function technique pro- 
vides a systematic way of understanding the transmission and conductance through short 
range scatterers in terms of resonant transport through double barrier structure. We par- 
ticularly point out to distinct regime of transport near the resonant value and off resonant 
values of such short range scatterrs. At and very close to the resonant value, the conductance 
in relatively large size sample shows very slow decay as a function of system size where as 
away from the resonant value the conductance in a large sample shows an algebraic decay 
as function of the system size with an exponent which is close to 0.5. Though our results 
are obtained by approximating the short range scatterers as delta function barrier, in order 
to obtain compact analytical expressions for transmittance and conductance, some of the 
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conclusions can be extended for more extended natured potential barriers as well. As our 
results suggested a transition from this resonant regime to the off resonant regime can be 
observed by introducing controlled disorder in graphene based superlattice structure and 
may well be used to suggest graphene based electronic devices. 
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